Size-Dependence of the Wavefunction of Self- Assembled Quantum Dots 
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The radiative and non-radiative decay rates of InAs quantum dots are measured by controlling the 
local density of optical states near an interface. From time-resolved measurements we extract the 
oscillator strength and the quantum efficiency and their dependence on emission energy. From our 
results and a theoretical model we determine the striking dependence of the overlap of the electron 
(_) ' and hole wavefunctions on the quantum dot size. We conclude that the optical quality is best for 

D , large quantum dots, which is important in order to optimally tailor quantum dot emitters for, e.g., 

^"*^ . quantum electrodynamics experiments. 
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Semiconductor quantum dots (QDs) have attracted 
significant attention recently as nano-scale light sources 
for all-solid-state quantum electrodynamics experiments 
P, 0, 0, 0, 0, 0] . Major advancements have culminated in 
the demonstration of strong coherent coupling between a 
single QD and the optical mode of a cavity 0,1. The 
coupling strength between the emitter and the cavity is 
determined by the oscillator strength, which is an intrin- 
sic property of the emitter. For atomic transitions the 
oscillator strength attains only discrete values depending 
on the choice of atom and is determined by the electro- 
static potential. In contrast, the QD oscillator strength 
can be ingeniously tailored due to the influence of size- 
confinement on the electron- hole wavefunction 0, 0, 0] • 
Consequently the oscillator strength can be continuously 
tuned by varying the size of the QDs. Surprisingly the ex- 
act size-dependence of the optical properties of the exci- 
ton has remained an open question. Understanding these 
excitonic optical properties is much required in order to 
optimally engineer QDs for enhanced light-matter inter- 
action. 

In this Letter we present measurements of the oscilla- 
tor strength of the ground-state exciton in self-assembled 
QDs. The detailed dependence on the QD size is mapped 
out by time-resolved measurement of spontaneous emis- 
sion at different emission energies. We employ the mod- 
ified local density of optical states (LDOS), caused by 
reflections in a substrate-air interface, to separate radia- 
tive and non-radiative decay contributions. This method 
was pioneered by Drexhage for dye molecules [l(| and 
used also to extract the quantum efficiency of Erbium 
ions 11 1 and colloidal nanocrystals [H, El . Here we use 
this method to accurately determine the dependence of 
the oscillator strength on the quantum dot size. The pre- 
cise measurements of the oscillator strength allows us to 
determine the size-dependence of the electron and hole 
wavefunction overlap. 



Time-resolved spontaneous emission is measured from 
a series of samples that contain identical ensembles of 
InAs QDs positioned at controllable distances to a GaAs- 
air interface. The wafer is grown by molecular beam 
epitaxy on a GaAs (100) substrate where 2.0 monolay- 
ers of InAs are deposited at 524° C followed by a 30 s 
growth interrupt and deposition of a 300 nm thick GaAs 
cap. The QD density is 250 um~ 2 . A 50 nm thick layer 
of AlAs is deposited 650 nm below the QDs for an op- 
tional epitaxial lift-off. The wafer is processed by stan- 
dard UV-lithography and wet chemical etching, whereby 
samples with different distances between the QDs and 
the interface are fabricated on the same wafer, see insert 
of Figure [TJ3. The distances z from the QD-layers to the 
interface are measured by a combination of secondary 
ion mass spectroscopy and surface profiling with typical 
precisions of ±3.0 nm. 

The QDs are excited by optical pumping of the wet- 
ting layer states at 1.45 eV using ~ 300 fs pulses from 
a mode-locked Tksapphire laser. The excitation spot 
has a diameter of ~ 250 u.m and the excitation density 
is kept at 7W/cm 2 with a repetition rate of 82 MHz. 
Under these conditions less than 0.1 excitons per QD 
are created, i.e., only light from the QD ground state 
is observed. The spontaneous emission is collected by 
a lens (NA=0.32), dispersed by a monochromator, and 
directed onto a silicon avalanche photo diode for time- 
correlated single-photon counting 14]. The detection en- 
ergy is varied between 1.17 eV and 1.27 eV to probe dif- 
ferent sub-ensembles of the inhomogenously broadened 
ground state. The spectral resolution of the monochro- 
mator is 2.6 meV, which is narrow relative to the band- 
width of the LDOS changes. The time-resolution of the 
setup is 48 ps given by the full width half maximum of the 
total instrument response function. All measurements 
are performed at 14 K. 

Figure [TJ\ shows the spontaneous emission decay for 



2 




100 150 200 

Distance [nm] 



FIG. 1: (color online). A Decay of the spontaneous emission 
recorded at 1.20 eV for two different distances to the inter- 
face of z = 109 nm (green, upper curve) and z = 170 nm 
(blue, lower curve). The bi-exponential fits (solid red lines) 
result in T f = 0.91ns" 1 , T s = 0.09 ns" 1 for z = 109 nm and 
r f = 1.15 ns" 1 , T s = 0.10 ns" 1 for z = 170 nm. The goodness- 
of-fit parameters Xr are respectively 1.17, and 1.11, close 
to the ideal value of unity [3| verifying the bi-exponential 
model. B Measured decay rates versus distance z to the 
GaAs-air interface (dots). Calculated LDOS projected onto 
a dipole orientation parallel to the interface (solid blue line). 
Calculated LDOS including dissipation at the surface (dashed 
red line). The inset is a schematic drawing of the sample. C 
Coefficient of correlation versus the number of data points 
excluded in the modeling of the data in B. 



QDs positioned at two different distances from the GaAs- 
air interface and recorded at an emission energy of 
1.20 eV. A clear change in the decay curve is observed 
with distance to the interface. The decay of the QD 
ground state is very well modeled as a bi-exponential de- 
cay, J(t) = Afe" r '* + A s e~ r,,t + C, over the complete 
time range of the measurement. The background level 
C is determined by the measured dark count rate and 
after-pulsing probability of the detector. The fast decay 
takes place on a time scale of about 1 ns corresponding 
to the decay of bright excitons in InAs QDs. The slow 
decay time is approximately 10 ns and does not system- 
atically depend on the distance and is most likely due to 
recombination of dark excitons [l5| . In the remainder of 
this Letter we will focus only on the fast decay rate. 

The decay rates measured at 1.20 eV are presented in 
Figure QJ3 as a function of distance from the QDs to 
the interface. A damped oscillation of the total decay 
rate with distance is observed. The data are compared 
to the LDOS calculated for GaAs (assuming n = 3.5) 



and projected onto a dipole orientation parallel to the 
interface (solid blue line). Only the parallel component 
is relevant since refraction in the interface reduces the 
solid angle for light collection, and additionally the QD 
orientation is predominantly parallel to the interface 16]. 
The measured decay rate T(w,z) is the sum of a non- 
radiative r nra d(w) and a radiative r ra d(w,,z) decay rate. 
The latter is proportional to the projected LDOS p(ui, z) 
and depends explicitly on the distance z to the interface 
and the optical frequency to. It is calculated as the sum 
over all available electromagnetic modes projected onto 
the orientation of the dipole, and is obtained using a 
Green's function approach [17| . We define the radiative 
decay rate for QDs in a homogeneous medium r|? a ° d n (cj) 
and express the measured total decay rate as 
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where Phom(w) is the LDOS of a homogeneous medium 
of GaAs. 

Excellent agreement between experiment and theory is 
observed in Figure [TJ3 for distances of z > 75 nm. This 
explicitly confirms the validity of the theoretical model 
used to extract prop erties of the emitter as opposed to 
a previous work [12| . For QDs closer than 75 nm to the 
GaAs-air interface the measured decay rates are system- 
atically larger than the calculated rates. We exclude that 
this effect is due to tunneling out of the QDs, which has 
been observed only within 15 nm from a surface (l8j . An 
increased non-radiative loss may be due to scattering or 
absorption at the surface of the etched samples. This 
dissipation in the surface is modeled as a thin absorbing 
surface layer, which creates an optical surface state. The 
dashed line in Fig. [TJ3 is obtained by including a 5 nm 
thick layer with refractive index of 3.5 + l.Oi, which leads 
to increased rates near the interface in agreement with 
the experimental data. 

For QDs sufficiently far away from the interface the 
influence of any surface effects is negligible, and our data 
can be used to reliably extract QD properties. We de- 
termine the data points that are not influenced by the 
dissipation in the surface as follows: Eq. fT]) reveals a 
linear relation between the measured rate and the calcu- 
lated normalized LDOS. We therefore perform a linear 
regression analysis and obtain the linear correlation pa- 
rameter \r\ as the data close to the interface are excluded 
point by point, cf. Figure [TJU. After excluding the seven 
closest data points the correlation parameter converges to 
unity, hence Eq. ([1]) is valid. By comparing experiment 
and theory we determine the radiative and non-radiative 
decay rates at 1.20 eV to be r^ a ° d n = 0.95 ± 0.03 ns" 1 
and r nra d = 0.11 ± 0.03 ns -1 , which is the most accurate 
result to date for QDs. 

The measurements have been performed for six differ- 
ent energies within the inhomogenously broadened emis- 
sion spectrum of the QDs. The inhomogeneous broaden- 
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ing reflects the different sizes of QDs such that small QDs 
correspond to high emission energies and vice versa. The 
radiative and non-radiative decay rates extracted from 
the measurements are plotted in Fig. The increased 
non-radiative recombination rate at higher energies could 
indicate that carriers can be trapped at the QD surface 
since the relative importance of the surface is large for 
small QDs. While such a size dependence would be gen- 
eral for all QDs, the absolute values of the non-radiative 
rates could depend on sample growth. Surprisingly the 
radiative rate is found to decrease with increasing en- 
ergy. This behavior is due to the decrease of the overlap 
between the electron and hole wavefunctions as the size 
of the QD is reduced, as discussed below. Our method 
allows to extract the size-dependence of the QD emis- 
sion without any implicit assumption of vanishing non- 
radiative recombination, as opposed to previous works 
8, 19]. In fact such an assumption would in our case lead 
to the incorrect conclusion of an increased wavefunction 
overlap with reduced QD size. 

Fermi's Golden Rule relates the radiative decay rate 
and the oscillator strength 
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where n is the refractive index of GaAs, u> is the fre- 
quency of the optical transition, m e is the electron mass, 
eo is the vacuum permittivity, q is the elementary charge, 
and Co is the speed of light in vacuum. For the QDs 
emitting at fkJ = 1.20 eV, the measured value of 
results in an oscillator strength of / osc = 13.0 ± 0.4. For 
comparison various estimates of the oscillator strength 
based on absorption measurements have been reported 
in the literature and are generally in the range of / OS c=5- 
10 [2^, 21 1 . However, the technique implemented here 



provides unprecedented precision since it only relies on 
accurate measurements of the distance of the QDs to the 
interface and is independent of, e.g., the QD density. Ad- 
ditionally, the quantum efficiency, i.e., the ratio of the 
radiative decay rate to the total decay rate, is extracted. 
We find QE = 90 ±4% at the emission energy of 1.20 eV, 
which confirms that a high quantum efficiency is feasible 
with an ensemble of emitters and not only with selected 
single QDs [12]. The intrinsically high QD quantum ef- 
ficiency can be further increased by tuning the size as 
discussed below. 

The energy dependence of the oscillator strength and 
the quantum efficiency are presented in Figure [3JA_. Both 
quantities are seen to decrease with increasing energy. 
The quantum efficiency decreases from around 95% to 
80% and the oscillator strength from 14.5 to 11 over the 
inhomogcnously broadened emission spectrum. This re- 
sult shows that large QDs with a high exciton confine- 
ment potential have much better optical properties than 
smaller QDs. Our results thus shed new light on the opti- 
mum design of solid-state QED experiments, and strong 
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FIG. 2: (color online). Left axis: Photofuminescence from 
the inhomogenously broadened ground state measured at z = 
281 nm (solid line). Right axis: Radiative (green triangles) 
and non-radiative (blue circles) decay rates versus emission 
energy. 
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FIG. 3: (color online). A Oscillator strength (triangles) 
and quantum efficiency (squares) versus energy. B Measured 
(squares) and calculated (red curve) energy dependence of the 
overlap of the envelope wavefunctions of electrons and holes. 



coupling was indeed observed for large QDs P, @, [1] • 

Our measurements provide new insight into the size 
dependence of the QD's wavefunctions. Thus, the spa- 
tial overlap between the electron and hole wavefunctions 
\^e,h{<jj)) can be obtained from the oscillator strength. 
Within the effective mass approximation, valid in the 
strong confinement limit when Coulomb effects are neg- 
ligible, the electron (or hole) wavefunction can be fac- 
torized in a conduction (or valence) band Bloch wave- 
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function for InAs | m c / v ) , and an electron (or hole) enve- 
lope function |F c /h(w))- The overlap of the electron and 
hole wavefunctions is related to the oscillator strength 



via 
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where e is the polarization unit vector of the electromag- 
netic field and p is the electron momentum. Evaluation 
of the matrix element is performed as an average over 
all possible orientations of the polarization vector. As 
strain lifts the degeneracy of the light-hole and heavy- 
hole bands, only transitions from the conduction band to 
the heavy-hole band are included. The result can be ex- 
pressed as |{it v |e-p|it c )| 2 = m P ,E n /6, where E n = 22.2 eV 
is the Kane energy for bulk InAs [23j. From Eq. ([3]) 
and the measured radiative decay rates, the wavefunc- 
tion overlap is obtained, see Fig. [3J3. The wavefunction 
overlap decreases from about 0.75 at 1.17 eV to 0.63 at 
1.27 eV. The reduction in the oscillator strength stems 
from the increased mismatch between the electron and 
hole wavefunctions with decreasing QD size. The reduc- 
tion of the wavefunction overlap is due to the more sen- 
sitive size-dependence of the electron wavefunction com- 
pared to the hole wavefunction due to their difference 
in effective mass. As the size of the QD is decreased, 
the electron will penetrate deeper into the barrier than 
the hole thus reducing the overlap [24j . The detailed un- 
derstanding of how the size affects the overlap between 
the electron and hole wavefunctions is crucial in order to 
optimally tailor QDs for efficient coupling to light. 

The observed energy dependence of the oscillator 
strength is compared to a simple effective-mass QD 
model. We use finite-element-method calculations to ob- 
tain the energy levels and the corresponding wavefunc- 
tions of the electron and holes. The curve in Fig. [3)3 
displays the wavefunction overlap calculated for a lens- 
shaped QD with a radius of 7 nm and a height varying be- 
tween 1.8 nm and 3.0 nm. The following realistic param- 
eters are used: a wetting layer thickness of 0.3 nm, 60% 
of the band-edge discontinuity is in the conduction band, 
and the GaAs content in the QDs is taken to be 25%. 
Good agreement with the experimental data is observed, 
and the theory clearly confirms a pronounced reduction 
of the electron-hole wavefunction overlap as the size of 
the QD is decreased. The general validity has been tested 
by calculating for different sizes, shapes, and amount of 
GaAs content in the QD, and they all show a decrease of 
overlap with increasing energy. The same behavior is also 
obtained from more involved QD models also in the pres- 
ence of strain 25j. Interestingly, the radiative rate was 



observed to increase with energy for colloidal nanocrys- 
tals in agreement with theory [19(, which is opposite to 
the results reported here for InAs QDs. This illustrates 
a striking difference in the optical properties of colloidal 
nanocrystals compared to self-assembled QDs, which is 



due to their different sizes and confinement potentials. 

In summary, we have measured the radiative and non- 
radiative decay rates of InAs QDs by employing the mod- 
ified LDOS near a dielectric interface. The oscillator 
strength and quantum efficiency of the QDs and their 
dependence on the emission energy were accurately de- 
termined. The radiative decay rate decreases with in- 
creasing energy leading to a reduction of the oscillator 
strength. In contrast the non-radiative recombination 
rate increases with increasing energy corresponding to a 
reduction of the quantum efficiency. Consequently QDs 
emitting on the low-energy side of the inhomogenously 
broadened ground state transition are most suitable as 
nanophotonic light sources due to their optimized optical 
properties. The experimental findings are explained by a 
model of the QD taking the size-dependence of the wave- 
functions into account. Our results demonstrate how QD 
wavefunctions can be tailored to achieve improved cou- 
pling to light, which is needed in order to take full advan- 
tage of the potential of quantum electrodynamics devices 
based on QDs. 
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